SuperWIMP dark matter and 125 GeV Higgs boson in the minimal GMSB 
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Recently, both the ATLAS and CMS experiments have observed an excess of events that could be the first ev- 
idence for a 125 GeV Higgs boson. We investigate an implication of the CP-even Higgs boson with mass around 
125 GeV in the context of the minimal gauge mediated supersymmetry breaking (mGMSB). In mGMSB, grav- 
itino is the lightest sparticle (LSP) and hence the dark matter candidate. We consider the so-called superWIMP 
scenario where the dark matter gravitino is non-thermally produced by the decay of the next-to-LSP (NLSP) 
bino-like neutralino after its freeze-out. For a given tan /3 and the number of the messengers (Nm ) fixed, we find 
that the rest of the mGMSB parameters, the SUSY breaking parameter and the messenger scale, are completely 
fixed by the conditions of simultaneously realizing the observed dark matter abundance and the 125 GeV Higgs 
boson mass, leading to the NLSP neutralino mass around 1.5 — 2 TeV and the gravitino mass around 3 — 7 GeV, 
depending on the values of tan /3 and Nm- The lifetime of the NLSP is found to be shorter than 1 sec, so that 
the success of the big bang nucleosynthesis remains intact. The non-thermally produced gravitino behaves as 
the warm dark matter with the free-streaming scale found to be Apg ~ 0.1 Mpc, whose value is reasonable for 
observations of the power spectrum on both large and sub-galactic scales in the Universe. 



The low energy supersymmetry (SUSY) is arguably one of 
the most promising candidates for new physics beyond the 
Standard Model (SM). The minimally supersymmetric exten- 
sion of the SM (MSSM) provides us with not only a solution 
to the gauge hierarchy problem but also various interesting 
phenomena such as the successful SM gauge coupling unifi- 
cation, the radiative electroweak symmetry breaking, predic- 
tion of the SM-like Higgs boson mass, and a candidate for the 
dark matter in our Universe. It has been expected that some 
of sparticles can be discovered in the near future, most likely 
at the Large Hadron ColUder (LHC). 

As well as the discovery of new particles, the discovery of 
the Higgs boson is another major goal of the physics program 
at the LHC, in order to confirm the origin of the electroweak 
symmetry breaking and the mechanism of particle mass gen- 
eration. Recently, the ATLAS [1] and CMS |2] experiments 
have reported an excess of evens that could be the first evi- 
dence of the Higgs boson with mass of around 125 GeV IH- 
The observations are supported by recent analysis of the Teva- 
tron experiments [3]. The 125 GeV Higgs boson has a great 
impact on SUSY phenomenology, because the MSSM has a 
prediction of the SM-like Higgs boson mass as a functions 
of soft SUSY breaking parameters, in particular, stop masses. 
Detailed studies for realizing the Higgs mass around 125 GeV 
is a current hot topic in SUSY phenomenology fT]. In most 
of the studies, the constrained MSSM (CMSSM) or sHght ex- 
tention of the CMSSM is adopted for the boundary conditions 
on the soft SUSY parameters at the scale of the grand unified 
theory (GUT). 

Although the CMSSM or more generally, the context of su- 
pergravity mediated SUSY breaking is a very simple bench- 
mark in examining the sparticle mass spectrum, this scenario 
potentially suffers from the SUSY flavor problem fl- The 
gauge mediated SUSY breaking (GMSB) |6] offers the com- 
pelling resolution for the SUSY flavor problem due to the 
SUSY breaking mediation via the flavor-blind SM gauge in- 
teractions. In this paper, we investigate an implication of the 
125 GeV Higgs in the context of the mGMSB. There were 



several studies on this context fl^, and very recently a more 
comprehensive studies with the parameter scan has been per- 
formed [8]. Although these studies identified a parameter 
space to realize the Higgs mass around 125 GeV, phenomenol- 
ogy of the gravitino dark matter has not been studied in detail. 
In this paper, we complete our study on the mGMSB by con- 
sidering not only the realization of 125 GeV Higgs mass but 
also the the gravitino dark matter phenomenology, in partic- 
ular, the so-called superWIMP scenario. We will show that 
once the values of tan /3 and the number of messengers {N,n) 
are fixed, the sparticle mass spectrum is completely deter- 
mined and the resultant sparticle mass spectrum is consistent 
with all phenomenological constraints. 

In the mGMSB, we introduce the superpotential for the 
messenger sector, 

where S* is a chiral superfield in the hidden sector with (S) = 
M + O^F, and $i and $i are a vector-like pair of messengers 
of the 5 + 5 representation under the SU(5) GUT gauge group. 
For simplicity, we use the SU(5) GUT notation throughout 
the paper. In order to maintain the successful SM gauge cou- 
pling unification, a pair of messengers should be in a complete 
SU(5) representation and have an SU(5) invariant mass term. 

Soft SUSY breaking terms can be extracted from the SUSY 
wave function renormalization coefficients with the thresh- 
old corrections by the messengers [9]. The MSSM gaugino 
masses at a scale /i < M is given by (we assume F ^ AP) 

M,(^) = ^AiV™, (2) 
47r 

where i = 1,2,3 coiTespond to the SM gauge interactions 
of SU(3), SU(2) and U(l)y, respectively, and A = F/M 
is the SUSY breaking parameter When we neglect Yukawa 
coupling contributions, the MSSM scalar squared masses at 



^ < M aTS given by 



mjifi) 



An 



A2iV„G,(/i,M), (3) 



whereG,(Ai,M) = ^f + ^{l-^f) with^ = a,(Af)/a,(M). 
Here bi are the beta function coefficients for different SM 
gauge groups, Ci are the quadratic Casimir, and the sum is 
taken corresponding to the representation of the sparticles un- 
der the SM gauge groups. The messenger has no contribution 
to 74-parameter at the messenger scale, A{M) ~ 0. In the fol- 
lowing numerical analysis, we employ the SOFTSUSY 3.3.1 
package |10] to solve the MSSM RGBs and produce mass 
spectrum, where the following free parameters are defined at 
the messenger scale: 



Nm, M, A, Cgrav, tan/3, sign(^). 



(4) 



For simplicity, we set /i > 0. In our analysis we assume 
Cgiav = 1, which means F is the dominant SUSY breaking 
source. In this case, gravitino mass is given by 



MA 



(5) 



where Mp = 2.4 x 10^^ GeV is the reduced Planck mass. 
Since flavor-dependent supergravity contributions to sfermion 
masses are estimated as Am^ ^ m%, we need to set m% ^ 

TO^, equivalently, M <^ 10^^ GeV, to make the flavor- 
independent GMSB contributions dominant. Therefore, grav- 
itino is always the LSP and hence the dark matter candidate. 
From Eqs. and ^ we see that the NLSP is most likely 
bino-like neutralino for N„i = 0(1). 

As is well known, stop loop corrections play a crucial role 
to push up the SM-like Higgs boson mass from its tree-level 
value, nifi ~ mz cos 2(3. Since the corrections logarithmi- 
cally depend on stop mass, a large stop mass of 0(10 TeV) 
is necessary to realize rrih — 125 GeV. This corresponds to 
A - 10^ GeV from Eq. © and mass of the NLSP bino-like 
neutralino to be m^j ^ 1 TeV from Eq. (|2|i, for Nm = 0(1). 

Now let us consider phenomenology for the gravitino dark 
matter This is quite different from the usual one with the dark 
matter as a weakly interacting massive particle (WIMP), be- 
cause gravitino couples to the MSSM particles super-weakly 
and (unless the gravitino is light, say mg, < 1 keV) it has 
never been in thermal equilibrium in history of the Universe. 
This fact brings uncertainty in evaluating the relic abundance 
of the dark matter and thus one may think this scenario less in- 
teresting. However, there is a very appealing possibility with 
a super-weakly interacting dark matter particle, the so-called 
superWIMP scenario flT, T?) (see also for the super- 
WIMP scenario in non-standard cosmology). In this scenario, 
the superWIMP dark matter is mainly produced via the decay 
of a long-lived WIMP, so that its relic abundance is given as 



V mx 



(6) 



where VI xh? would be the thermal relic abundance of the 
long-lived WIMP (X) if it were stable. In this scenario, nev- 
ertheless the superWIMP has never been in thermal equilib- 
rium, its relic abundance is calculable as in the usual WIMP 
dark matter scenario. 

We adopt the superWIMP scenario in the mGMSB, where 
the superWIMP is the LSP gravitino and the WIMP X is 
the bino-like neutralino. Note that bino-like neutralino, if it 
were the LSP, tends to be over-abundant because of its small 
annihilation cross section. Therefore, the suppression factor 
jTig/m^j ^ 1 can work well to adjust VloKih^ to be the ob- 
served value lil4]. With a fixed tan j3 and Nm, we first cal- 
culate particle mass spectrum by SOFTSUSY 3.3.1 package 
for various values of M and A. Then, the relic abundance of 
the NLSP neutralino is calculated by using the micrOMEGAs 
2.4.5 J3] with the output of SOFTSUSY in the SLHA for- 
mat [16]. Multiplying by the factor 111^/771^, we finally ob- 
tain the relic abundance of the gravitino dark matter. 
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FIG. 1 : Contours on the (A/, A)-plane for tan /3 = 10 and A^™ = 1. 
The dark shaded region (in red) and the horizontal shaded region (in 
yellow) satisfy the constraints of the observed dark matter abundance 
(0.1064 < n^h^ < 0.1176) and the SM-like Higgs boson mass 
(124 GeV <mh < 126 GeV), respectively. 



Fig. 1 shows our numerical results in the (Af , A)-plane, for 
tan /? = 10 and Nm = 1. In the dark shaded region (in red), 
the resultant gravitino relic abundance is consistent with the 
observed value (vf\, n^K^ = 0.1120 ± 0.0056. The parame- 
ter sets in the horizontal shaded region (in yellow) predict the 
SM-like Higgs boson mass in the range of 124 GeV < rah < 
126 GeV. The dashed lines correspond to the gravitino mass 
ttlq — 1, 2, 3, 4, and 5 GeV, respectively, from left to right. 
The parameter region which simultaneously satisfies the ob- 
served Higgs boson mass and the dark matter relic abundance 
can be pined down by the overlap of the two shaded regions, 
(A/, A) ~ (9.2 X 10^2 GeV, 1.2 x 10*^ GeV). 

Here we give a semi-analytical explanation of our results. 
The relic abundance of the bino-like neutralino is well approx- 
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imated as 118 1 



.7 X 10-" GeV"^ (n + l)x^ 



n+l 



(7) 



where n 



1, Xf 



20 is the freeze-out tem- 



perature, and cr^ is the pair annihilation cross section of the 
bino-like neutralino through the exchange of the right-handed 
charged sleptons in the ^-channel, 



3g'^ tan* 0yyr{l + r^) 
27rm?^(l+r4) 



(8) 



with r = (Afi/mg^)^. Using Eqs. (|2|i and (O, we can check 
that the factor r(l + r^)/(l + r*) is almost constant 0.12 
for various values of Nm = 0(1) and M = lO^^-i* QeV. 
Thus, we find 



Ml 



«m^^^(xMA2. (9) 



Since sfermion masses logarithmically depend on the messen- 
ger scale (see Eq. (|3]l), the predicted Higgs mass determined 
by stop masses is almost independent of M — 10^^^^* GeV. 
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FIG. 2: Contours on the (M, A)-plane for tan /3 = 10 and Nm = 
1,2,3, 4. The solid lines correspond to the region satisfying Q^h^ — 
0.112, from left to right for Nm ~ 1, 2, 3, 4, respectively. Along the 
dashed lines the SM-like Higgs boson mass is predicted to be rrih ~ 
125 GeV, from top to bottom for A'^^ ~ 1, 2, 3, 4, respectively. 



We repeat the similar numerical analysis for various val- 
ues of Nm and tan (3. Fig. 2 depicts the results for = 
1, 2, 3, 4, with tan/3 ~ 10. The solid lines show the contours 
along which fl^^h^ — 0.112 is satisfied, from left to right cor- 
responding to Nm = 1,2,3,4, respectively. All lines are very 
close to each other except that the line for Nm = 4 is substan- 
tially deviating from the other lines for A/ < 5 x 10^"^ GeV. 
This is because for such a parameter region, the NLSP bino- 
like neutralino becomes more degenerate with lighter stau as 
M is lowered, so that the co-annihilation process of the NLSP 
with the stau becomes more effective and reduces the abun- 
dance. In order to compensate this reduction and achieve the 
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TABLE I: Particle mass spectra (in GeV) for various Nm and tan /?. 
The fundamental parameters in the mGMSB, M and A, have been 
fixed so as to satisfy two independent conditions, Q,Qh^ = 0.112 
and TUh ~ 125 GeV. 



observed relic abundance, a relatively higher neutralino mass, 
equivalently a larger A is needed. The dashed horizontal lines 
correspond to the Higgs boson mass prediction niH — 125 
GeV, from top to bottom for Nm — 1, 2, 3, 4, respectively. 
We find that stau becomes the NLSP for Nm > 5. The anni- 
hilation process of the NLSP stau is very efficient, so that the 
suppression factor m^/nix makes the gravitino abundance 
too small. Hence, we do not consider the case with Nm > 5. 

In Fig. 2, we can pin down the parameter sets {AI, A) which 
simultaneously satisfy the two conditions, fl^h^ = 0.112 and 
nih = 125 GeV, at the intersections of the solid and dashed 
lines, for each Nm value. In Table 1, we list particle mass 
spectra for Nm = 1,2,4. Here we also list the result for 
Nm = 1 and tan (3 ~ 45 as another sample. 

In the superWlMP scenario, the NLSP has a long lifetime 
and its late-time decay is potentially dangerous for the success 
of big bang nucleosynthesis (BBN). In our scenario, the NLSP 
is bino-like neutralino and its lifetime is estimated as 1. 1 1.1 



0.74 sec X 



( y 

VI GcV/ 



1 TeV 



(10) 



for its main decay mode B j + G. For the parameter 
sets we have determined, the formula leads to < 1 sec. 
Therefore, the NLSP neutralinos decay before the BBN era 
^ 1 sec and the success of BBN remains intact. 
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The gravitino dark matter is produced by the decay of the 
bino-like neutrahno and carries a large kinetic energy when 
it is produced. Such a non-thermally produced dark mat- 
ter behaves as a "warm" dark matter, rather than the cold 
dark matter (CDM). While the CDM scenario has made great 
success in the structure formation of the Universe at large 
scales > 1 Mpc, the recent high resolution N-body simula- 
tions showed that the CDM scenario predicts too much power 
on sub-galactic scales to be consistent with the observations. 
It has been pointed out [l^] that this problem can be solved 
by a non-thermally produced dark matter if its comoving 
free-streaming scale (Aps) at the time of the matter-radiation 
equahty (tEQ — 2.0 x 10^^ sec) is around 0.1 Mpc. 

The free-streaming scale of the gravitino dark matter can be 
calculated as lEoll 



A 



FS 



dt ~ 2woiEQ(l + zeq) 



log 



^'o(1 + ^eq)^ wo(1-I-2;eq) 



,(11) 



where zeq — 3000 is the red shift at ^eq, and vq ~ 
is the current velocity of the gravitino dark matter with the 
present temperature of the cosmic microwave background Tq, 
the temperature of the Universe Tj ~ y/1 scc/t^ x lO^'^ 
Ge V, and the energy of the gravitino Ej = 'ti^ /2 when the 
gravitino is produced from the NLSP neutralino decay. It is 

easy to find ug — 1.0 x 10~^ ( I^IeX j ^ which is indepen- 
dent of the gravitino mass, and the free-streaming scale is 



Aps ^ 0.18 Mpc X 



1 TcV 



3/2 



(12) 



The mass spectra listed on Table 1 give the free-streaming 
scale Aps — 0.1 Mpc, which is nothing but the value suit- 
able for solving the problem with the CDM scenario on sub- 
galactic scales. 
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